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MULTI-USER ADAPTIVE ARRAY RECEIVER, AND METHOD 

DESCRIPTION 

5 

TECHNICAL FIELD: 

The invention relates to a receiver system comprising an antenna and a receiver, the 
antenna comprising an array of antenna elements. It also relates to the receiver per se and 
to a receiving method. The invention is especially, but not exclusively, applicable to array 
10 receivers and methods for use in base stations of digital cellular telecommunications networks 
or access points of wireless local area networks (LANs). 

BACKGROUND ART: 

The invention is applicable in systems wherein multiple users simultaneously ihake use 

15 of a common carrier or use distinct carriers with some amount of overlap/crosstalk between 
them such as (i) carrier reuse-within-cell (RWC), also called Space-Division Multiple Access 
(SDMA) because of the need for an antenna array to spatially discriminate co-channel signals; 
(ii) Code-Division Multiple Access (CDMA) systems where multiple users transmit in the 
same band using distinct codes; and (iii) Time-Division Multiple Access (TDMA) and/or 

20 Frequency-Division Multiple Access (FDMA) systems where users are not perfectly separable 
in time and/or frequency, i.e. they interfere with one another either in time (e.g. because of 
dispersive channels) and/or in frequency (e.g. because of excess bandwidth due to imperfect 
channel filtering) thus leading to adjacent-channel interference (ACI). 

Mathematical expressions in this patent specification are based upon complex 

25 baseband notation. 

Array antenna radio receivers typically are employed at the base stations or access 
points of digital communications systems (e.g. mobile telephone networks, broadband 
wireless access for Internet and/or wide-area networking, etc.) to improve reception link 
quality (i.e. provide robustness against multipath fading) and/or reduce interference levels, 

30 where interference can include thermal noise and man-made signals which exist in the desired 
signal's band. Since such systems typically accommodate large numbers of simultaneously 
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active users in any given cell or sector, the base station receiver must be capable of 
maintaining a plurality of radio links. 

Known antenna array radio receiver systems comprise an array of antenna elements 
coupled to a signal receiving section (also referred to as a radio-frequency (RF) front-end) 
5 which in turn is coupled to a signal processing section. The signal receiving section processes 
the branch signals from the different antenna elements independently, in separate branches, 
and performs on each branch signal standard downconversion, demodulation, filtering to 
isolate the channel of interest and, possibly, some transformation on the signal to bring it to 
a form usable by the signal processing section (e.g. analog-to-digital conversion if the signal 
1 0 processor is digital) . The signal processor takes the information from all of the branches (i.e. 
the demodulated, filtered and suitably transformed signal data from each individual antenna 
element) and, using one of a number of appropriate known techniques, combines and 
processes it to extract a useful signal y(t), which is the best possible estimate of the desired 
user signal s 0 (f). 

15 In the context of wireless communications, the received vector x(t) (i.e. the received 

signal across all array elements) is made up of a desired signal s Q (f) transmitted by a "desired 
user's" wireless terminal, interfering signals s m (f) transmitted by competing terminals which 
operate in the same frequency band or in adjacent bands with some amount of crosstalk being 
present, and white noise n(t). Hence, in non^dispersive (i.e. narrowband) environments 

20 *(/) = c 0 (f> 0 (/) + £ c m ®sj® + *{t). (1) 

m=l 

where €^(0 is an Nxl vector of complex elements describing the channel from the mth 
terminal to all of theN array elements, Mis the number of interfering signals, n(t) is the white 
thermal noise vector, and c 0 (t) is an Nx 1 complex vector describing the channel from the 0th 

25 terminal which, by convention, is that of the desired user. 

In such a context, the iunction of the antenna array radio receiver is to isolate the 
desired user signal s 0 (t) from the interferers and white noise as well as compensate for 
distortions introduced in the channel c 0 (i) (e.g. multipath fading) so that, at all times, the 
array output signal y(t) approximates the desired user signal s 0 (f) as closely as possible. 

30 Typically, the receiver combines the branch signals from the individual antenna 

elements simply by means of a linear weight-and-sum operation. If an ^-element array is 
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considered andx(/) is the Nxl vector of the array element outputs, the array output is defined 
as 

y(t) = w*(ftx(& (2) 
where w(f) is the Nxl complex weight vector and (-) H denotes the Hermitian transpose (i.e. 
5 complex conjugate transpose of its argument, be it a vector (as it is in the above) or a matrix). 

Although it is time-varying, the weight vector varies slowly compared to the input and 
output signals, since it tracks changes in the channels, not in the signals themselves. When 
a combiner operates according to equation (2), it is termed a linear combiner and the entire 
receiver is designated a linear array receiver. 
10 Typically, the receiver collects statistics of the input signal x(f) and uses them to 

derive a weight vector which minimizes some error measure between the array output y(t) 
and the desired signal s 0 (f) One of the most common error measures in such applications (i.e. 
adaptive filtering) is the mean-square error 

* - (CKO -*o«] 2 ) = ([^(0^(0 -*o(0P) (3) 

15 which forms an ^-dimensional quadratic surface with respect to the weight vector elements. 
The minimization of this criterion forms the basis of minimum mean-square error (MMSE) 
linear array receivers (also called optimum combiners). 

(Note: Henceforth, the dependence upon time / in equations will be omitted for the 
sake of clarity.) 

20 Adaptive filtering theory indicates that the best combination of weights in the MMSE 

for a given sequence of received data is 

" = «2*o ( 4 ) 
where R**is the covariance matrix of the received array outputs and is given by 

* (5) 
25 where (•> denotes the expectation (i.e. the ensemble average) of its argument. 

Such array receivers are suitable for use where time dispersion due to multipath 
propagation does not extend significantly beyond a single symbol period. That is, there is 
little or no intersymbol interference (ISI). 

When the channels carrying useful signals do exhibit significant ISI, the traditional 
30 solution is to use an equalizer, which is an adaptive filter whose purpose is to invert the 
channel impulse response (thus untangling the ISI) so that the overall impulse response at its 
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output will tend to be much shorter in time and have an ideal, flat (or equalized) frequency 
spectrum. 

The signal processing portion of the standard linear equalizer works in the same way 
as a linear adaptive array receiver except that the signal sources, i.e., the elements of the input 
5 vector*, are not points in space (i.e. the array of antenna elements) but points in time. The 
signals are tapped at a series of points along a symbol-spaced delay line (termed a 
tapped-delay line or TDL), then weighted and combined 

While the implementation of the signal processing apparatus forboththe equalizer and 
the array receiver can be identical (minimization of the MSE by adaptive weighting of the 

10 inputs), the performance will differ. Because signals are physically sampled at different points 
in space by the array receiver, it is very effective at nulling unwanted signal sources or 
co-channel interference (CCI). However, it has limited ability against intersymbol interference 
(ISI) due to dispersive, i.e. frequency*selective , fading, since the latter is spread in time. On 
the other hand, the equalizer is adept at combatting ISI but has limited ability against CCI. 

15 In environments where both ISI and CCI are present, array reception and equalization 

may be combined to form a space-time array receiver. The most general form of the latter 
is obtained when each weighting multiplier in a narrowband array receiver is replaced by a 
full equalizer for a total of N equalizers. Again the implementation of the signal processing 
section will be similar and will rely on equation (2) supra. The only difference is that the 

20 weight vector w and the input vector jc will each be longer. Indeed, for an equalizer length 
of L taps and an array size of N elements, the vectors w and x will both have LN elements. 

The canonical linear mean-square-error minimizing space-time receiver (i.e. the most 
obvious and immediate linear space-time receiver structure and also in certain respects the 
most complex) comprises an antenna array where each array element output is piped to a 

25 finite impulse response (FIR) adaptive filter, which in this context is referred to as an 
equalizer. Each adaptive filter comprises a tapped-delay line where taps are spaced by a 
symbol period or a fraction of a symbol period. For good performance, the length of the 
tapped-delay line should be equal or superior to the average channel memory length. In many 
cases, the number of taps this implies can be very large (e.g. 10-100 per adaptive filter). 

30 The weights multiplying each tap output must be constantly adapted to follow the 

changes in the channel(s) characteristics. This can be performed in various ways, either with 
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continuous or block-based adaptation and with or without the support of known training 
symbols. In most known systems, the weights are computed on a block-by-block basis (block 
adaptation) and each block contains a sequence of known training symbols for that purpose. 
In digital wireless communications, the block used for adaptation purposes will typically 
5 correspond to a data packet as defined by the networking protocol in use. Moreover, the 
channels can be considered static over the length of a block (i.e., the length of a block is 
significantly smaller than the channel correlation time). 

By adapting the weights to minimize a global performance index, e.g. the mean-square 
error between the desired signal and the S-T receiver output, the receiver usually performs 
10 the following: 

reduces or eliminates inter symbol interference (ISI) caused by frequency- selective 
fading in wideband channels; 

reduces or eliminates co-channel interference (CCI) from nearest cells where carriers 
are reused or from inside the cell (since the space-time processor permits reuse of carriers 
15 within cell - or sector - thanks to its power of spatial discrimination - often referred to as 
space division multiple-access (SDMA)); 

• improves output SNR (due to the array's larger effective aperture). 

Since wireless systems are typically interference-limited (i.e., interference is the main 
impediment which prevents capacity increase - accommodating more active users - above a 
20 certain limit), the first two benefits of space-time processors are mainly of interest in order 
to increase capacity 

To achieve maximal benefit, it is better to combine the S-T array with carrier 
reuse-within-cell (RWC). A number of previous patents disclose arrays (see , for example US 
5,515,378 and US 5,592,490) or space-time systems (see, for example, US5,828,658) 

25 applied in an SDMA (i.e. RWC) context. In such a system, separate S-T processors will 
have to be implemented for every user (all processors share the same physical antenna array 
and front-end receiver circuitry but have distinct equalizers and combiners). However, the 
base station has information (received symbols, channel characteristics) available about in-cell 
interferers since each in-cell interferer is another local S-T processor's desired signal. 

30 S-T processor architectures can be formulated to exploit this multiuser information 
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by establishing some type of connectivity between individual S-T processors to achieve one 
of two goals: 

• improve performance (reduced bit-error rate, improved interference nulling, etc.); 

• reduce complexity and cost. 

5 It is known to exploit multiuser information to perform "joint detection" of many 

users, for example by constructing a global multiuser MSE criterion, thus improving 
performance of an array receiver (with respect to single user detection) at the cost of 
increased complexity [2], [3]. 

It is also known that, with appropriate space-time processing, it is possible to combine 

10 SDMA with adequate temporal processing to mitigate the intersymbol interference (1SI) 
present in wideband dispersive channels [6]. 

One of the main disadvantages of previously-known space-time processing receivers 
is their great complexity and cost, especially if multiuser detection is employed and/or 
temporal processing employed. 

15 It is known to reduce bandwidth requirements in forward-channel probing transmitters 

by tracking only long-term variations in the channels (i.e., the subspace structure) [1 1] but 
that approach is not applicable in receivers without seriously limiting user capacity. 

DISCLOSURE OF INVENTION: 

20 The present invention seeks to at least mitigate the disadvantages of such known array 

receiver systems and to this end provides a multiuser space-time array receiver, and system 
incorporating same, exploiting multiuser information in order to reduce complexity and cost. 

According to one aspect of the present invention, there is provided an array receiver 
for processing signals received from a plurality of transmitting users via an array antenna 

25 having an array of N antenna elements providing a set of antenna signals (x iy x 2 >— , x N ), 
respectively, each comprising information from each user, characterized by a common 
preprocessing section for sampling each of the antenna element signals (x ly x^,.-., x N ) and 
processing the samples of at least some of said antenna element signals to form a plurality of 
subspace basis signals, for example subspace signals, (y 0 > . . j>m) together having fewer space- 

30 time dimensions than the space-time dimensions of the combined antenna signals, and a 
plurality of signal processing units each having a plurality of inputs coupled to the common 
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preprocessing unit for receiving all of the basis signals, each processing unit processing and 
combining said basis signals to produce a respective one of a set of estimated received signals 
( 2 o> ■ > *m) eac h f° r a corresponding one of the users, the common preprocessing section 
comprising filtering means for combining all of the antenna signals (x x > x j7 ..., x N ) to provide 
5 said plurality of basis signals (y^. y^, each of the basis signals comprising a different 
combination of the antenna signals, each of the signal processing units combining the basis 
signals to provide a user-specific output signal, and updating means for periodically updating 
parameters of the filtering means used for deriving each particular basis signal such that each 
user-specific output signal will exhibit a desired optimized concentration of energy of that 

10 user's received signal as received by the array antenna. 

In preferred embodiments, the updating means comprises means for adjusting said 
parameters in dependence upon channel characteristics of all user channels. Each of the 
processor units then may comprise means for weighting the basis signals before combining 
same, the weights being adjusted in dependence upon channel characteristics of all user 

15 channels, wherein the parameters of the filtering means are updated less frequently than the 
weights of the processors. 

The number of basis signals may be equal to the number of desired users. 
The updating means may comprise a training sequence generator for generating a 
training sequence for the corresponding user, co variance matrix estimation means responsive 

20 to the training sequence and the antenna signals for providing a covariance matrix embodying 
long-term statistics for the channel of that user, and eigenvector estimation means for 
extracting from said covariance matrix at least the dominant eigenvector, elements of said 
dominant eigenvector being applied to said filtering means as weights for updating said 
parameters. 

25 Preferred embodiments of the first aspect of the present invention address the 

complexity issue by (1) creating, from the antenna element outputs, a common basis of filters 
(i.e. filter bank) useful for all users 1 processors, (2) adapting this basis based on the 
slowly-varying statistical channel structure, thus reducing the computational burden, and (3) 
by selecting for each user only a few (e.g. 2 or 3) most significant filter outputs from the 

30 common basis for rapid adaptation. 
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Preferably, when the array receiver system is employed in a space-division multiple 
access (SDMA) communications system, the plurality ofbasis filters in the preprocessing unit 
and the plurality of subsequent receiver signal processing units are both assigned to an 
ensemble of transmitting users sharing a common channel (Le. frequency band) at the same 
5 time. 

Alternatively, the plurality ofbasis filters and subsequent receiver sections could be 
assigned to an ensemble of transmitting antennas belonging to the same user, yet transmitting 
different bit sequences in order to provide a higher aggregate bit rate. This latter 
configuration corresponds to a multi-input (MIMO) link. It should be understood that, in the 

10 following, references to a "user 75 in a SDMA context would translate to 'transmitting 
antenna" in a MIMO context and that the techniques described in a SDMA context otherwise 
are directly applicable in a MIMO context. 

According to a second aspect of the invention, there is provided a receiver for 
receiving signals from a plurality of transmitting users via an array antenna having an array 

15 of N antenna elements providing a set of antenna signals (x u x 2t x^ respectively, each 
comprising information from each user, said receiver characterized by a common 
preprocessing section followed by a plurality of receiver sections, each corresponding to a 
different one of the users and coupled to the outputs of the common preprocessing section, 
the preprocessing section sampling each of the antenna signals (x x , * 2 ,..., x N ) and processing 

20 the samples of at least some of said antenna element signals to form a plurality of basis 
signals (y c ,---. Jm) together having fewer space-time dimensions than the space-time 
dimensions of the combined antenna signals, and a plurality of signal processing units each 
having a plurality of inputs coupled to the common preprocessing unit for receiving all of the 
basis signals, each processing unit processing and combining said basis signals to produce a 

25 respective one of a set of estimated received signals (z 0 ,.,., each for a corresponding 
desired one of the users, 
the common preprocessing section comprising 

(i) means for maintaining through periodic updates a set of dominant subspace filters, 
each of which being matched to one of the users of interest, and the outputs of which being 
30 used by the subsequent receiver sections, to be processed and combined in order to yield an 
estimate of the desired signal for each user of interest; 
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(ii) means for periodically estimating and/or updating the component weights of the 
dominant subspace filters by correlation, with a known training sequence or with the user's 
spreading code in a CDMA system or with any other signal strongly correlated with the user 
of interest's signal, in combination with appropriate temporal averaging to isolate 

5 subspace-level information, as opposed to instantaneous channel characteristics, and 

(iii) means for periodically or dynamically estimating and/or updating the component 
weights (and/or any other parameters of interest) of the receiver sections fed from the 
preprocessing section in a manner and at a rate such that instantaneous channel changes are 
tracked to provide a reliable and consistent estimate of the desired signal. 

10 Preferably, in embodiments of either aspect, when the array receiver is employed in 

a code-division multiple access (CDMA) communications system, the plurality of basis filters 
forming the common basis and the plurality of subsequent receiver signal processing units are 
both matched to: 

(1) an ensemble of users sharing the same spreading code, if code re-use is employed in 
15 the said communications system, or; 

(2) an ensemble of users with different codes, in which case the array receiver system can 
further separate the users 1 signals and possibly compensate discrimination problems due to 
code correlation, power capture, etc. 

In a CDMA system, the usual despreading can be performed at the outputs of the 
20 basis filters. Alternatively, the spreading operation at the transmitter can be considered a 
part of the radio channel, in which case it is natural to despread in the basis filters. 

The array receiver of either aspect could also be employed at the base station of a 
time-division multiple access (TDMA) communications system or a frequency-division 
multiple access (FDMA) communications system which does not employ carrier re-use. In 
25 such a case, the plurality of dominant basis filters of the common preprocessing unit and the 
plurality of subsequent receiver signal processing units are both matched to an ensemble of 
users which are not perfectly separable in time and/or frequency, i.e. they interfere with one 
another either in time (e.g. because of dispersive channels) and/or in frequency (e.g. because 
of excess bandwidth due to imperfect channel filtering) thus leading to adjacent-channel 
30 interference (AO). 
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Other aspects of the invention include an array receiver system comprising a receiver 
according to the first or second aspect, in combination with a said array antenna, and the 
corresponding method of operating the array receiver. 

Thus, according to a third aspect of the invention, there is provided a method of 
5 receiving signals from a plurality of transmitting users via an array antenna comprising an 
array of AT antenna elements providing a set of antenna signals (x l3 x N ), respectively, 
each comprising information from each user, the method comprising the steps of receiving 
signals from a plurality of transmitting users via an array antenna having N antenna elements 
providing a set of antenna signals (x 1? x 2> ..., x w ), respectively, each comprising information 
10 from each user, and being characterized by the steps of: 

sampling each of the antenna signals; 

preprocessing the samples of at least some of said antenna element signals (x u x 2 , •> 
x^) to form a plurality of basis signals (y 0 ,,. M y^ together having fewer space-time 
dimensions than the space-time dimensions of the combined antenna signals, and 
15 processing and combining said basis signals (y 0 , . . Jm) t0 produce a set of estimated 

received signals (2 0 v> *m) eac h for a corresponding one of the users, 
the preprocessing including the step of 

combining all of the antenna signals (x 1? x^..., x N ) to provide said plurality of 
basis signals (y 0 , . . y^) such that each of the basis signals comprises a different combination 

20 of the antenna signals, 

the processing and combining step comprising the step of combining the basis signals 

(y 0 , to provide a series of user-specific output signals, 

the method further comprising the step of periodically updating parameters used for 
deriving each particular basis signal such that each user-specific output signal will 
25 exhibit a desired optimum concentration of energy of the received signal of that 

particular user as received by the array antenna. 

Preferred embodiments of this third aspect of the invention comprise method steps 
corresponding to the functions of embodiments of the array receiver of the first and second 

■ 

aspect. 

30 In one preferred embodiment, the receiver: 
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(1) maintains, through periodic updates, a set of dominant basis filters, each of 
which is matched to one desired user among the users of interest, and the outputs of which 
are processed by the subsequent receiver sections and combined in order to yield an estimate 
of the desired signal for each desired user; 

5 (2) periodically estimates and/or updates the component weights of the dominant 

subspace filters by correlation, with a known training sequence or with the user's spreading 
code in a CDMA system or with any other signal strongly correlated with the user of 
interest's signal, in combination with appropriate temporal averaging to isolate subspace-level 
information, as opposed to instantaneous channel characteristics; 

10 (3) periodically or dynamically estimates and/or updating the component weights 

(and/or any other parameters of interest) of the receiver sections fed from the prefiltering 
section in a manner and at a rate such that instantaneous channel changes are tracked to 
provide a reliable and consistent estimate of the desired signal. 

In embodiments of any of the first, second and third aspects of the invention, the 

15 receiver may comprise a series of standard linear MMSE space-time processors (i.e. one 
possible embodiment of the receiver sections), one for each of the M+l signals, operating on 
the transformed input vector y[n]. However, this method only results in a net reduction of 
numerical effort if the number of signals M+l taken into account is significantly lower than 
the number of antenna elements N (in which case the dimensionality of the input vector is 

20 reduced from N to MM). 

Embodiments of any of the three aspects of the invention may include space-time 
matched filtering. This provides a much greater potential complexity reduction and makes 
the invention more widely applicable. Thus, to further decrease computational cost, a logical 
extrapolation of the above concept is to extend the eigenfiltering to the temporal - as well as 

25 the spatial - domain. In this case, only M+l taps are left to be actively adapted (at every 
packet) for each user (as opposed to NL taps for a conventional system where N is the 
number of elements and L is the required adaptive filter length. To achieve acceptable 
performance, it is normally required that MzN\ therefore this system will reduce the number 
of actively adapted taps by at least a factor of L.) 

30 It so happens that a large portion of the I SI will in most cases be handled by the 

eigenfiltering. Indeed, the angle spread of an impinging signal at the base station is typically 
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narrow due to the height of the base, i.e. most scattering activity then occurs in the 
immediate vicinity of the subscriber. This has the effect of making the covariance matrix of 
the signal under consideration heavily biased towards the first few eigenvalues [1 1 J. A bank 
of primary eigenfilters eliminates the ISI associated with the first eigenvalue. Furthermore, 
5 it has been shown that a memoryless combiner (such as those that follow the eigenfilter bank) 
has some ability to reduce ISI [13]. 

In cases where these two ISI reduction steps are not enough to warrant satisfactory 
performance, more dimensions can be added to the dominant subspace space-time filters to 
eliminate the ISI and CCI associated with the secondary and subsequent eigenmodes at the 
10 cost of increased complexity, since more taps will have to be actively adapted in the receiver 
sections. 

According to another embodiment of the invention, the receiver preprocessing 
sections can include adaptive equalization, thus reducing the need in certain cases for a large 
number of subspace dimensions to adequately handle the ISI. 
15 The foregoing and other objects, features, aspects and advantages of the present 

invention will become more apparent from the following detailed description, taken in 
conjunction with the accompanying drawings, of preferred embodiments of the invention, 
which are described by way of example only. 



20 BRIEF DESCRIPTION OF THE DRAWINGS: 

Figure 1 is a simplified block schematic diagram of an array receiver system having 
a receiver and an array antenna comprising an array of antenna elements, 

Figure 2 is a more-detailed block schematic diagram of a part of the receiver showing, 
in more detail, a dominant subspace filter for one user; 

25 Figure 3 is a flowchart depicting computation of updated weights for the dominant 

subspace filter; 

Figure 4 is a flowchart depicting computation of principal eigenvector estimates for 
use in updating the weights; 

Figure 5 is a flowchart depicting computation of secondary eigenvector estimates; 
30 Figure 6 is a block schematic diagram of a modified receiver system; and 
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Figure 7 is a flowchart depicting adaptation of weights to changing channel 
conditions. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS : 
5 To facilitate understanding of the construction and operation of the preferred 

embodiments, some basic theory will first be presented. 

As discussed, the optimal MMSE solution can be obtained through a linear 
combination of all signals' matched-filters [7], [8]. Given an antenna array and a dispersive 
(i.e. ISI-inducing) propagation environment, it follows that the optimal MMSE solution can 

10 be obtained as a linear combination of all signals' space-time matched-filters. Here, a 
space-time filter matched to a given signal is a bank of N temporal filters, each of which 
processes one of the N antenna elements 1 outputs, and whose outputs are combined to 
maximize the said signal's power with respect to white noise and disregarding interference 
from the other man-made signals. 

15 This is advantageous in a multi-user SDMA context since the set of matched filters 

form a common basis which can be reused to obtain each signal's MMSE solution. In 
standard optimal architectures , independent combiners (sets of weights) must be maintained 
for each user and they must typically be recomputed from scratch at the start of a new packet 
because of the changing interference patterns. If a way can be found to maintain with low 

20 computational cost a matched filter for every active connection, computing an MMSE 
solution for a given user and packet becomes simply a matter of selecting the appropriate 
matched filters (corresponding to the active interferers in the packet of interest) and using 
their outputs as inputs to standard MMSE processors (one per desired user) which are 
adapted using the training sequence prefixes. The complexity of this approach is appealing 

25 when the system is designed in such a way that the number of inputs to the standard MMSE 
processors is substantially reduced with respect to a system in which the antenna elements 1 
outputs are directly processed. 

One method to approximate the behaviour of a matched-filter without having to track 
the multipath fading is to identify dominant subspaces of the users' vector channels. The said 

30 subspaces will contain most of the useful information yet vary at a much slower rate than the 
channels themselves. A dominant subspace is the reduced-rank space spanned by the few 
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most significant dimensions of the long-terra eigenstructure of the channel of interest. One 
case of interest (because it minimizes complexity) consists in dominant subspaces with a 
single dimension. When the corresponding eigenvector is used as a filter, the resulting device 
is termed an eigenfilter. 
5 To obtain better estimates of the desired signals, it may be necessary to perform 

eigenfiltering with subspaces having more than one dimension. The required number n of 
dimensions for a given level of performance is a function of the propagation environment. 
Eigenfiltering over multidimensional subspaces is termed hereafter dominant subspace 
filtering. 

10 Referring to Figure 1, an array antenna receiver system for receiving signals from a 

plurality of user transmitters (not shown) comprises an array antenna having a plurality of 
antenna elements, specifically TV elements 22/1, 22/N, each coupled to a respective one of 
a corresponding plurality of RF front-end processing units 26/1, 26/N of an RF receiver 
section 26, which units treat the signals from the antenna elements to produce N signals x i9 

15 . x N , respectively. Each of the RF front-end units 26/1, .. 26/N has its output coupled to 
the input of each of a set of A/+1 filters specifically subspace filters, 40/0, 40/M of a 
common preprocessing section 40. Each of the subspace filters 40/0, 40/M is matched 
to a respective one of an ensemble of A/+1 transmitting users, and has its output coupled to 
the input of each of a corresponding plurality of user-specific signal processors 60/0, . , . , 60/M 

20 of a signal processing section 60. Each of the signal processors 60/0, . . , 60/M processes the 
respective set of the subspace signals y 0 , y M of the subspace filters 40/0, 40/M, 
respectively, to produce a corresponding one of a plurality of estimates z 0> z M of the 
signals s 0? s M transmitted by the A/+1 users. 

The RF "front-end" units 26/1, .... 26/N are identical and of conventional 

25 construction, so only one will be described, with reference to the inset diagram of Figure 1 , 
As shown inset in Figure 1, RF front-end unit 26/N comprises a low-noise amplifier (LNA) 
28/N, a RF to IF downconverter 30/N, a channel filter 32/N (which isolates the required 
channel and rejects out-of-band noise and interference), and an analog-to-digital converter 
34/N for performing bandpass sampling. Alternatively, the IF or RF signal could be down 

30 converted to baseband prior to analog-digital conversion. The various alternatives and 
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compromises possible here are known to practitioners of the art. The output of the A/D 
converter unit 34/N is also the output of the RF front-end unit 26/N and is coupled to the 
input of each of the dominant subspace filters 40/0, 40/M. 

In each of the following embodiments, all of the dominant subspace filters 40/0, 
5 40/M are identical; although their structure differs from one embodiment to another. 

In the embodiment of Figure 1, the dominant subspace filters 40/0, 40/M are 
principal space-time eigenfilters. Since they are all identical, only the generic structure of the 
filter for a user m will now be described, with reference to Figure 2. 

Although the performance analysis will be presented in the frequency domain, the 
10 actual implementation can be made in the time domain. The eigenfilters then take the form 
of banks of ATtapped-delay lines 50/m Xj ...50//w N each with a series of one-symbol delays, the 
number of such delays being chosen to give a delay line length according to the typical 
memory length of channels in the band of operation. In each tapped delay line, a series of 
multipliers extract the delayed signals from respective taps of the delay line and multiply each 
15 of them by a respective complex weight. For example, in delay line 50m 1? having individual 
delays 52m u , ... 52m lL , a series of multipliers 54m u ,... 54m 1L multiply the tapped signals by 
complex weights w Ut w iu respectively, while, in delay line 50m N having individual delays 
52m Nl ,„. 52m NL , a series of multipliers 54m m ,.. 54^ multiply the tapped signals by complex 
weights w N1 , ... tvhl, respectively. The other tapped delay lines are similar. 
20 The outputs of the delay lines 50/rn^..., 50/m N , i.e., the signals from the multipliers 

54m 1L ,.,., 54^ respectively, are combined by a summer 52/m to form>> wb the primary 
eigenfilter output for user m. It should be noted that there can be any number of such 
eigenfilters whose combined outputs will make up the dominant subspace filter output i.e. 
subspace signal y, n . Thus, y m = \y m%l , . . y m J where \i is the number of eigenfilters defining the 
25 dimensions of the dominant subspace. This estimate y m is supplied to all of the signal 
processors 60/0,..., 60/M (Figure 1). 

The weights are dynamically adapted according to second-order statistics as will now 
be described. A training sequence generator 42/m generates a replica of user m'$ known 
training sequence in synchronism with reception of said training sequence as part of a 
30 received packet. The output s m [k] of the training sequence generator 42/m is used by a 
covariance matrix estimator 44/m to estimate user m's current channel covariance matrix Rj 
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across all space-time tap positions. This is then used to update a running estimate of user 
m'$ long-term covariance matrix S M . An eigenvector estimator 46/m then uses the estimate 
of H m as a basis to estimate its principal eigenvector. The individual components of said 
eigenvector constitute the current complex weights w ll7 w ih3 ... which are supplied 
5 to multipliers 54 n> 54 lL , ... 54^ in filter bank 48/m for use in subsequent processing of 
the received signals. 

The adaptation procedure used by the eigenfilters 50/^, ...50//7i N will now be 
described in detail with respect to the flowchart in Figure 3, where the joint operation of a 
bank of 8 eigenfilters (corresponding to 8 transmitting users) is detailed. 

10 Adaptation of the eigenfilters requires that a running estimate of each signal's 

long-term covariance matrix be maintained. This estimate could be updated every time a 
packet containing known training symbols is received from the user of interest. Since the 
long-term statistics change relatively slowly, however, the estimate update frequency for a 
given user is going to be much lower than the frequency of occurrence of the training 

15 sequence (which typically is provided in every packet from the user of interest). 

In step 3.1, the dominant subspace filter 40/m waits for the current estimation interval 
to elapse (where "estimation interval" refers to the relatively long interval for long-term 
estimation as discussed above) and v in step 3.2, waits for the start of the next time slot. 
Assuming that the known training sequence is a prefix and is thus at the start of said time slot, 

20 in step 3 .3 the portion of the received signal corresponding to the training prefix is stored in 
a buffer for further processing. 

Given that £ m is the "delay-extended" NLxl vector representing the space-time 
signature of user m over all eigenfilter taps, i.e. 



25 where djn} = SJfiT) is a sample of the vector impulse response of user m*s channel at the 
array input taken at delay «r(wth multiple of the symbol period), in step 3,4, the dominant 




(6) 



subspace filter 40/m obtains an estimate^ of c 0 the vector impulse response for user 0 



30 
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where s 0 [k] is user O's training sequence (obtained from the training sequence generator 42/0 
which is K symbols long and x[k,ri] is the space-time received signal vector over all 
eigenfilter taps corresponding to the fah training symbol of the mh training sequence. For 
example, for a fixed estimation interval of T 9 seconds and a symbol period of T seconds, 
5 = x(nT c + kT). (8) 

In step 3.4, the eigenfilter covariance matrix estimator 44/m computes an estimate of 
the covariance matrix for user 0 for the current interval according to 

3>M - jzi,WJL*r. (?) 

10 

By definition, the long-terra delay-extended covariance matrix for signal (s m ) is 

- (io) 

15 Accordingly, a running estimate of £ 0 (for user 0) is updated in step 3.5. The 

estimator 44/m obtains a running estimate of % m according to the following recursive 
relation. 

£>] = Y^-l]^(y-l)^ w [/i] s (11) 

20 where S m [»] is the NLxNL long-term covariance matrix estimate after processing of the mh 
received training sequence for the signal from user m, s m [k] is the kth symbol in user m's 
training sequence and y is a forgetting factor chosen as a function of training update 
frequency and the rate of change of the long-term covariance matrix in the propagation 
environment of interest. It is likely that y would take on a value between 0.8 and 1 in most 

25 systems. 

The covariance matrix estimator 44/m supplies the estimate ll D to eigenvector 
estimator 46/m which uses it in step 3.6 to estimate the principal eigenvector, thereby 
completing the weight estimation/update procedure for filter bank 48/0. 

In this embodiment, the estimation of the principal eigenvector is performed using the 
30 iterative power method [12]. This requires an initial estimate of the eigenvector and an 
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5 



estimate of the covariance matrix (obtained in steps 3.5, 3.10, 3.14). Based on these 
estimates, the dominant eigenvector of can be obtained according to: 

*2m = tMrtFw. («) 



where A„[n] is the estimate in the nth estimation interval of the dominant eigenvector over 
all NL S-T eigenfilter taps after the rth iteration of the power method (including 
normalization, i.e. |<2«[n]| 2 te 1). The convergence rate of the power method depends on 
the ratio i^2y where k x and X z are the largest and second largest eigenvalues of 2 W [»]- For 

10 a well-conditioned matrix and any arbitrary starting vector tf£\ri\ 9 convergence will normally 
occur within 10 iterations. 

Upon network entry of a new user, a large number of iterations (50-100) might be 
necessary to guarantee a good estimate of the dominant eigenvector regardless of the 
eigenvalue distribution. Afterwards, however, since the successive covariance matrix 

15 estimates SJ/f] vary little from one to the next, only a few (perhaps even 1 or 2) iterations 
of the power method will be required between covariance matrix updates. 

The eigenvector estimation procedure (performed by the eigenvector estimators 46/0, 
46/7) will now be described with reference to the flowchart in Figure 4 for user m. 
In step 4. 1 , the estimator 46/m compares estimation interval index n with 0; if n = 0, 

20 then the first estimation is being performed since the group of users of interest has entered 
the network. Therefore, there is no previous estimate of the user m's principal eigenvector 
A m0 [n] and an arbitrary estimate is used in step 4.2 to set the initial starting point tf^oM- 
In step 4,3, the number of iterations / is set relatively high (50) since the starting point is not 
necessarily close to the real eigenvector. 

25 If ri>0 in step 4.1, then, in step 4.4, the estimator 46/m sets the starting point A^\n\ 

to the eigenvector estimate A^ 0 [n- 1] obtained in the previous estimation interval. In step 
4,5, it sets the number of iterations / to 5. 

In step 4.6, the estimator sets the iteration index I to 0 and then, in step 4 .7, performs 
a first iteration of the power method according to 

30 = tMftnl (13) 
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In step 4.8, the estimator normalizes vector^toututytoyiddarefined estimate \n\. 
In step 4.9, the estimator verifies whether the prescribed number of iterations has been 
performed. If not, in step 4. 10, it increments the iteration index / and repeats steps 4.7 to 
4.10 as indicated by the loop back to step 4.7. This process repeats itself until / = 7-1. 
5 Finally, the final product of the iterative process A^~ l \n]k assigned as the current principal 
eigenvector estimate A m [n] . 

Every time a new update of the dominant eigenvector is obtained, its components are 
immediately transferred into the weight registers of the eigenfilter banks 50/m l3 50/m„ 
(Figure 2). 

10 The same procedure is repeated in steps 5,8 through 5.14 for users 1 to 7, i.e. filter 

banks 48/1 through 48/7. 

In an alternative embodiment, corresponding to that shown in Figure 3 but with the 

addition of the items shown as dotted lines and boxes, secondary eigenvectors are also 

computed in steps 3 .7, 3 . 1 1 , . . 3 . 1 5, adding a second output to all dominant subspace filters 
15 40/0, . 40/M and thus providing more flexibility to the subsequent signal processors 60/0, 
60/M This can provide better performance against intersymbol interference (I SI) and 

against co-channel interference (CCI) and/or lessen the requirement for temporal processing 

in processors 60/0, 60/M, as explained below. In fact, any desired number R of 

eigenvectors can thus be computed to achieve the desired cost/performance compromise 
20 and/or the desired complexity balance between the common preprocessing section 40 and the 

per-user processors 60/0, 60/M. 

Estimation of secondary and further eigenvectors must be performed in order of 

decreasing eigenvector importance. After estimating the principal eigenvector according to 

Figure 4, the estimator subtracts its contribution from the covariance matrix < The 

25 resulting covariance matrix, designated^, has a principal eigenvector which is approximately 

equal to the secondary eigenvector of £ W M; the latter can therefore be estimated by 

following the procedure described with reference to Figure 4. 

This procedure is detailed in Figure 5 for an arbitrary number R of dominant 

eigenvectors. In step 5. 1, the eigenvector order index r is set to 1 7 indicating the principal 
30 eigenvector. In step 5.2, the initial covariance matrix^ is set to . An estimate of the 

principal eigenvector u m [ri] is then obtained in step 5.3 based on A r and according to the 
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procedure outlined in Figure 4. The eigenvector estimator 46/m then verifies in step 5.4 
whether all R eigenvectors have been computed. If not, in step 5,5, it computes an estimate 
of the rth order eigenvalue according to 

5 

and, in step 5,6, subtracts the corresponding eigenvector fromy4 r according to 

A„ x -Ar-tf^tfr (15) 

10 It then increments the index r in step 5.7 and repeats the steps 5.3 to 5.6 for the subsequent 
eigenvectors. 

The output signals y } ,...,y M from the dominant subspace filters 40/0, 40/M, i.e., 
the outputs of the common preprocessing section 40, are used by per-user signal processors 
60/0, 60/M to provide user-specific estimated received signals z 0v .,, 2 M , respectively, 
15 which are estimates of the A/+1 desired signals, each signal processor using the outputs of 
all of the dominant subspace filters 40/0, 40/M to produce its respective estimate. 

It will be appreciated that the common preprocessing section 40 described with 
reference to Figures 2 to 5 yields a set of signals with a reduced number of dimensions for 
further processing. This set of signals, or basis, is adapted through long-term adaptation since 
20 it tracks only the subspace structure of the channels, not their instantaneous behaviour. 

This long-term loop (which corresponds to Figure 3) need only be performed once 
every tenth of a second. This estimation interval can correspond to several hundred packets. 
On the other hand, the method associated with the per-user processors 60/0,... 60/M is a 
short-term loop (to be described hereafter for the preferred embodiment with reference to 
25 Figure 7) which typically must be performed once per packet. 

The signal processors 60/0,... 60/M can take a number of forms. According to this 
preferred embodiment where the dominant subspace filters 40/0, 40/M are space-time 
principal eigenfilters, as described previously, the signal processors 60/0, 60/M simply 
comprise weight-and-sum structures across the eigenfilter outputs, as illustrated in Figure 6. 
30 Referring to Figure 6, signal processors 60/0, 60/M are identical so only one will 

be described, namely signal processor 60/0. It comprises a plurality of multinliers 64/0^ 
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64/0 M which apply weights w^, . . . , *>om to subspase signals y^ . . . , y M , respectively, from the 
dominant subspace filters 40/0, 40/M. As will be described later, the weights w^, 
are derived in dependence upon substantially instantaneous channel characteristics, and 
updated. Means for deriving and updating these weights would be known to persons skilled 
5 in this art and so, for purposes of clarity, this means is not shown in Figure 6 

The weighted signals are then summed by combiner 62/0 whose output z 0 is fed to the 
detector 80/0 (see also Figure 1). It should be noted that each of the other signal processors 
40/1, 40/M also uses all of the output signals > 0 , y M to obtain its respective one of 
outputs z u z M . 

10 It is also important to note that, should transmission be momentarily interrupted (such 

as in bursty data communications), no problem occurs since the weights are estimated afresh 
in every interval. 

In step 7 J, the eigenvector estimator 46/m (Figure 2) waits for the next packet to 
start. Upon acquiring a packet, in step 7.2 the estimator stores the portion of the packet 
15 which corresponds to the training sequence in a buffer (not shown). In step 7.3, it computes 
the estimate of the (M+l)x(M+l) short-term covariance matrix Ryy of the eigenfilter bank's 
output vector^ as follows: 



20 



K ■ |f>™*]" (16) 



No training sequence is necessary since, at this point, there is no need to discriminate between 
individual transmitted signals. Furthermore, the same matrix is used by all users so this 
step need only be performed once per loop if all users 1 training sequences are synchronized. 
In step 7.4, the estimator sets user index m to 0 and in step 7.5 it estimates user w's 
25 (Af+l)xl signature vector over the eigenfilter bank's outputs using the known training 
sequence transmitted by user m, i.e. 



30 The training sequence will be provided by the training sequence generator 42/m which 

is part of eigenfilter 40/m (Figure 2). The generator 42/m will then occasionally (when there 
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is an intersection with the longer estimation interval) feed the "long-term" estimation in 
dominant subspace filter 40/m and the "short-term" estimation in processor 60/m 
simultaneously. 

Indeed, although the "short-term" loop of Figure 7 typically will be performed on the 
5 order of 100 times more often than the "long-term" loop of Figure 3, both can use the same 
training sequences, presumably available in every packet. However, the long-term loop will 
obviously use fewer of them. 

It should be noted that the "short-term" covariance matrix and the channel estimates 
all should preferably be estimated over the same set of received samples (which correspond 
10 to training sequences sent simultaneously by all users) to ensure statistical consistency and 
prevent serious performance degradation. 

In step 7.6, the estimator computes the weight vector w M which minimizes the 
mean-square error according to 

= (is) 



15 



processor 



In step 7.7 the estimator 46/m transfers the weights obtained thereby to the signal 
;or 60/m. The estimator 46/m then determines, in step 7.8, whether the weights in all 



* * K 7 ■ ^ CP — — — 

signal processors have been updated. If they have, the estimator returns to step 7. 1 to await 
the next time slot and repeat the procedure. If they have not, the estimator increments index 
20 m in step 7.9 and repeats steps 7.5 through 7.8 to update the next processor. Steps 7.5 to 
7,9 are repeated until all of the processors 60/0, . . 60/M have had their weights updated and 
step 7.8 finds m=M, 

The embodiment described above with reference to Figure 6 is designated "space-time 
eigenfiltering followed by MMSE combining" . 
25 To recapitulate, embodiments of the present invention can be described in algorithmic 

fashion as comprising two loops, vis. long-term and short-term. The long-term loop can be 
summarized as follows. 
For every user m, 

(i) The short-term covariance matrix of user m l s signature over all NL taps of eigenfilter 
30 m is estimated on the basis of a known training sequence transmitted by user m\ 
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(ii) The short-term estimate is used to update a running estimate of the 
long-term-averaged covariance matrix of user ro's channel estimate (eqn. 1 1); 

(iii) Using the running estimate of the long-term covariance matrix and the estimate of its 
dominant eigenvector from the previous iteration as a starting point, said running estimate 

5 is updated by performing one or more iterations of the power method (eqn. 12). If secondary 
eigenfilters are implemented, the same procedure applies for updating the secondary 
eigenfilter except that the dominant eigenvector is a priori subtracted from the covariance 
matrix estimate. 

(iv) The computed weights (i.e. elements of the estimated eigenvector(s)) are transferred 
10 to the /nth dominant subspace filter. 

(v) The start of the next long-term training interval is awaited, and steps (i) to (iv) then 
are repeated. 

The short-term loop can be summarized as follows: 
1 . The (Af+ 1 )x(A/+ 1 ) short-term covariance matrix JJ^ of the eigenfilter bank's outputs 
15 is estimated; 
For every user m, 

2, User m x $ (M+l)xl signature vector is estimated over the eigenfilter bank's 
outputs; 

3. The weight vector w m = A'^i m is computed; 
20 4. The weights are transferred to the /wth combiner. 

5 . The start of the next short-term training interval (next packet transmission by 
mth user) is awaited, and then the loop (steps 2 to 4) is repeated. 

Alternative embodiments 

25 The signal processors 60/0, , . „ 60/M could also include equalization, thus performing 

space-time processing. Such an extension is relatively straightforward to one skilled in the 
art and has the advantage of improving the performance (at the cost of additional complexity) 
in terms of signal quality and/or alleviating the need for many subspace dimensions in the 
preprocessing section 40 in order to obtain a given level of performance, 

30 In an alternative embodiment, the eigenfilter banks 48/0, 48/M perform strictly 

spatial processing, leaving all temporal processing to the per-user signal processors 60/0, 
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60/M. The eigenfilter banks 48/0, 48/M then have the structure depicted in Figure 2 but 
with only one column of weights, i.e. L=l. Such an embodiment is designated "spatial 
eigenfiltering followed by MMSE space-time processing". 

It should be noted that the latter embodiment provides better performance gains but 

5 a smaller complexity reduction than "space-time eigenfiltering followed by MMSE 
combining". In fact, a net complexity reduction compared with respect to conventional 
space-time processing can only be obtained if the number of antenna elements is greater than 
the number of users. However, this is usually the case even in conventional space-time 
processors in order to provide some gain against multipath fading in addition to spatial 

10 discrimination of users 1 signals. 

While it is general practice to assume that the channels can be considered static over 
the length of a block (i.e., the length of a block is significantly smaller than the channel 
correlation time), the present invention is applicable equally well in other cases where 
continuous tracking (using adaptive algorithms such as the least-mean-square (LMS) or the 

15 Kalman filtering algorithm) is necessary. 

If, in fact, continuous tracking is implemented, it may not be necessary to provide 
frequent training sequences. Indeed, both subspace filtering and weight computation updates 
can be performed using past decisions as training symbols, provided the latter are reliable 
("decision-directed adaptation"). Training sequences, while less frequent, would still be 

20 required to: (1) initialise the system when a new link was formed so that its first decisions 
would be reliable enough to start the tracking procedure; and (2) periodically reset the system 
to minimize errors due to tracking. 

Blind adaptation techniques could also be used, in which case training sequences 
would not be required at all. Likewise, the principles of the invention apply equally well to 

25 analog waveforms as opposed to digitally-modulated signals. 

The transmitting stations need not be limited to using a single antenna. If they have 
multiple antennas, thus forming multiple-input, multiple-output (MIMO) links, embodiments 
of the invention as described here can be modified appropriately in a number of ways while 
retaining the essence and advantages of the invention. For example, each transmitter antenna 

30 element belonging to the same user could have at the receiver its own dominant subspace 
filter. Thus, an ensemble of dominant subspace filters would feed a single per-user signal 
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processor which could perform standard MIMO reception techniques such as layered 
space-time (LST) successive cancellation. 

Error-correction coding, whether unidimensional or bidimensional (in MIMO links), 
can also be incorporated in ways that should be obvious to a skilled practitioner of the art. 
5 Likewise, a variety of alternatives to linear MMSE processing can be considered for 

the per-user processors 60/0 7 60/M without departing from the scope of the invention. 
Possibilities include decision-feedback processing, delayed decision-feedback, multi-user or 
MIMO decision-feedback, maximum-likelihood sequence estimation (MLSE), etc. 

The basis signals matched to each user could be formed using alternative techniques 
10 which are not based on the subspace structure of the channels. They could, for example, be 
based on estimates of the main directions-of-arrival characterizing each user's signal. 

The invention can also be applied to CDMA systems. Thus, for example, the usual 
despreading could be performed at the outputs of the subspace filters 40/0... 40/M. 
Alternatively, a bank of despreaders could be provided at the input of the preprocessing 
15 section 40 and supply all of the despread signals to each of the dominant subspace filters. 

Complexity reduction 

For the purpose of comparing complexity, an example will be considered of a lOMb/s 
system with packets of 68 bytes (roughly the size of an ATM ceil). A guard byte is inserted 

20 between each pair of successive packets. If there are 8 users (i.e. M-7) who send packets 
simultaneously once every ten slots on the same carrier, since there are 181 15.94 slots per 
second, the users of interest are transmitting at a rate of 1 8 1 1 . 59 packets per second. At this 
rate, channels typically will be sufficiently different from one packet to the next due to 
multipath fading to warrant retraining of the per-user processors 60/0. . . 60/M at every packet, 

25 It will also be assumed that all adaptive filters have a length of 10 symbol-spaced taps; each 
packet contains a known training sequence of 32 bits; and the array has 10 elements (i.e. 

The long-term covariance matrix is assumed to have a worst-case 90% correlation 
time of 0.5 s [9]; its estimate will be estimated every 0.1 s and the power iteration will also 
30 be performed every 0.1 $. Furthermore, it is assumed that diagonal loading (i.e. adding a 
small constant to all elements of the diagonal) is used whenever it is necessary to invert a 
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covariance matrix larger than the length of the training sequence (i.e. larger than 32x32) 
since it might otherwise be singular under these conditions 

Table 1 illustrates the relative numerical complexity of the two proposed structures 
versus standard MMSE space-time processing in terms of the number of multiplication and 



5 addition operations required. AJI figures are totals for all 8 users. 







standard MMSE 


spatial eigenfilter + 
S-T MMSE 


S-T eigenfilter + 
MMSE comb. 


10 


long-term 
adaptation 
(per iteration) 


8267200 mult. 
8260550 adds. 


15320 mult. 
14800 adds. 


1413200 mult. 
1372000 adds. 


15 


short-term 
adaptation 
(per iteration) 


1. 498- 10 10 mult. 
1. 496- 10 10 adds. 


4271360 mult 
4266840 adds. 


7808 mult. 
7644 adds. 




TOTAL 
(per second) 




7.73 8- 10 9 mult 
7.730- 10* adds. 


2.828- 10 7 mult. 
2.757- 10 7 adds. 



20 Table 1 : Relative numerical complexity of proposed structures compared with conventional 
MMSE space-time processing. 



THEORY OF OPERATION 

While not wishing to be limited by theory, an explanation of the theoiy of operation 
25 will now be given to facilitate understanding of the preferred embodiments. 

Given an eigendecomposition of the long-term average covariance matrix of the signal 
(i.e. its subspace structure), the eigenvector corresponding to the largest eigenvalue 
constitutes what will be called here the primary eigenfilter. Formally, the NxN long-term 
correlation matrix (where N is the number of receiving antenna elements at the base station) 
30 of signal $ m (t) (transmitted by the /wth out of M+ 1 users) in a flat-fading environment can be 
defined as 

E„ = {xjto'jtfl (19) 

where 

*„(') = C J<)SM (20) 

35 

and xjt) is the received signal from user m, c^r) is the N x 1 baseband equivalent vector 
channel between user m and the array (also called user tn's spatial signature), sjt) is the 
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useful transmitted signal from user m and the expectation <■> can be interpreted as either the 
time average over a period of time long enough to eliminate short-term channel fluctuations 
or an ensemble average over the distribution of possible channel realizations. 

It is well-known that, in almost all terrestrial propagation environments narrowband 
5 (i.e. flat fading) wireless channels can be accurately represented in the short-term as either 
zero mean (Rayleigh-type fading) or non-zero mean (Rician-type fading) complex gaussian 
variables. It follows that the vector c m taken at any time instant is a complex gaussian vector 
characterized by its long-term covariance matrix (which is equal to H m in theRayleigh-fading 
case) and its mean vector u where 

is the general definition of the long-term covariance matrix of user /w's vector channel- 
Without loss of generality, the remainder of this description will assume that the 
fading is Rayleigh and user m l s covariance matrix can thus be denoted by S w without 
15 ambiguity. 

For frequency-selective fading channels, the correlation matrix can be defined as a 
frequency-dependent matrix: 

£ w GO - ] rieJ&Xtr), (22) 

20 

where &\\%\ denotes the Fourier transform taken over the delay r variable and <•>, denotes 
averaging over the time variable /. Also, the dispersive channel impulse response c(t,-z) is 
defined as the echo received at time t + x originating from an impulse sent at time t. 

Consideration will be given first to the primary spatial eigenfilter in the general case 
25 of frequency-selective fading channels . For spatial filtering only, development proceeds from 
a "frequency-flat" covariance matrix for signal m, obtained by further frequency- or 
delay-averaging: 

2 m7 -4t- r~E* co4T- <* m w m o,*r> w (23) 

J\ max J mm 
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It should be noted that although the analysis presented here is in the frequency 
domain, the preferred embodiments described hereinafter are implemented in the time domain. 

The two definitions for £ m (in the delay domain and the frequency domain) can be 
proven equivalent through Parseval's relation. 
5 The primary spatial eigenfilter for signal m is then simply the eigenvector of S m 

corresponding to its largest eigenvalue. Its length is equal to the number of antenna elements 
and it is implementable as a set of weights used in combining the outputs of the array. 

Statistically, it can be shown to be the fixed combination of weights providing the 
highest average signal output without tracking the multipath fading. Some analytical and 
10 simulation results have indicated that the long-term correlation matrix changes relatively 
slowly even with mobile subscribers and can in general be assumed fixed for periods of the 
order of a second [9], This assumption has also been exploited to form the basis of downlink 
beamforming systems in [10] and [11]. In the broadband wireless context, it is reasonable to 
expect that the rate of change would be even slower since the subscribers are fixed. This 
1 5 implies that, in all cases, the eigenfilters can be computed in the background u sing a long-term 
tracking adaptation system demanding a negligible numerical effort This is where the 
complexity advantage of this invention lies, as described in the "preferred embodiments" 
section. 

If there are M+l signals and the desired signal is signal Sq, let us define a spatial 
20 transformation on the input vector which can take the form 

y[n] Ux[h] = £/S[w] **[/!], (24) 
where 6[n] is the Kronecker delta function, x[n] is the array input vector at time index n and 



25 U = 



H 
H 

h 



H 



(25) 



is an (M+l)xN matrix with u m being the primary spatial eigenfilter of the mth signal as 
described above. At the output of this transformation, the correlation matrix of the mth signal 
30 is expressed 

K V)-UY, m <f) U H for all fe (26) 
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and the desired signal signature at the output of transformation U becomes 



(27) 



where c 0 (f) is the desired signal signature prior to the transformation. 
5 It follows that the performance of an MMSE space-time combiner using the vector 

y[n] as input can be analyzed using known techniques and formulas for standard linear 
MMSE adaptation parametrized on the modified covariance matrices. 

The rest of this section comprises a performance analysis in the frequency domain of 
the embodiment with space-time subspace filtering, without adaptive equalization and further 
10 assuming without loss of generality that the dominant subspace filters have one dimension. 

Development proceeds by defining the prefiltering transformation (i.e. the eigenfilter 
bank) corresponding to the primary eigenvectors on "frequency-extended" space-time 



15 



25 



30 



vectors: 



(28) 



where 



£ = 



min 



f. 



*nth 




fin/re 



/ J 



(29 



is theNN b xl frequency-extended array input vector obtained by splitting the band of interest 
20 into N b bins of width A b significantly smaller than the coherence bandwidth (i.e. the fading 
can be considered flat within a single bin). Likewise jns the (M+^NbXl frequency- extended 
vector of the eigenfilter outputs and is expressed 



= A, 



/. 



/. 



min 



A. 

/ — - 



\1 



/ J 



(30) 



where each element of vector p is a complex number representing a superposition of flat 
fading channels and y(f) is the M+lxl frequency-dependent output vector of the space-time 
eigenfilter bank. Also, it should be noted that y is vector of size (M+l)N b xl . 
Also, 



(31) 
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where is the frequency-extended form of the principal space-time eigenfiJter 
corresponding to signal s„. 

All those quantities are time-varying (except for U which is fixed in the short-term) 
but time dependence is omitted for clarity. 

The system performance can be assessed by defining "frequency-crunched 0 channel 



vectors: 



m 



Lin* 



(2k- 1)\ 



l 6> 



(32) 



10 where A fc = ""^ ls ^ e width of each aforementioned frequency bin and djf) is an 

(M+l)xl vector defining user m's signature after the transformation t/'such that its frequency 

extended version = U8 . Furthermore, £ is a frequency-extended vector defined in the 

same manner as in eqn. (29). 

The above operation can be rewritten as a further linear transformation of the form 

!5 = VK (33) 

where V is an (M+l)x(M+l)N b matrix of the form 



20 



V- 



V I 

1 0 0 (A/-1 zeros) ... 10 0 (M~l zeros) 

0 10 (M-l zeros) ... 0 10 (Af-1 tferos) 

0 0 1 (A/^l reras) .„ 0 0 1 (M-l zeros) 



■ ■ 



(34) 



and d m is the frequency-extended version of the signature d m The corresponding short-term 
covariance matrix is given by 

25 K»*£*JZ+*u.i<> 2 (35) 



The output of a linear MMSE combiner taking y as input but which does not perform 
temporal processing (i.e. equalization) can be written 

z = (36) 
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Since only a single weight multiplies the output of each space- time filter, the extended 
weight vector is spectrally flat. In other words, 

w(f) = w. (37) 
Therefore the extended weight (M+l)N a xl vector has the form 
5 tiP = [^y,...,^]. (38) 

This implies that the output can be expressed 

z = w H y, (39) 



10 



where M , 

r A 



\ 



(40) 



This can also be written as a linear transformation 

y = ^ (4i) 

From MMSE filtering theory, the optimal weight vector in this context is 

"opt = *M*o (42) 
1 5 Therefore, the optimum MMSE performance of space-time eigenfiltering followed by a linear 
MMSE combiner is 

Ho = _ (43) 

By virtue of the transformation in equation (33), it is obvious that d 0 is a complex 
gaussian vector with covariance matrix 

20 ^ S 0 = VUJj q V H V h 9 (44) 

and that is equivalent to a short-term interference-plus-noise covariance matrix obtained 
from a fictional M+l element array in a flat fading environment. Each of its interference 
terms is obtained from a complex gaussian vector parametrized on appropriately transformed 
covariance matrices as per equation (44). 

25 While this approach reduces complexity considerably compared to standard 

space-time processing by exploiting a multiuser channel basis, it does so at the cost of 
reduced performance. Specifically, the robustness against fading (ISI) will be reduced 
somewhat. In cases where this is not acceptable, the system can be augmented by the 
addition of a second dimension to the subspaces (thus implementating secondary eigenfilters) 

30 or more dimensions to the subspace filters. 

The output of the secondary bank of eigenfilters is obtained, as in equation (28), by 
defining a prefiltering transformation corresponding to the secondary eigenvectors 
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(45) 



where 



and is the principal space-time eigenfilter 



U 2l 



1M 



(46) 



corresponding to signal (s^). 



It follows that in a system extended to include a secondary eigenfilter bank, the 
overall output vector of the prefiltering section is defined: 



10 



$ - 



9 m 








■ 


„ ?i 






■ 

* 


■ ■ 















X. 



(47) 



15 In this context, user m's signature becomes 



Ml 



ft 



(48) 

after the prefiltering transformation. Likewise, user m's frequency-crunched signature 
20 becomes: 

d = , (49) 




The rest of the development is identical to the case where only the primary eigenfilters 
are used. Similar generalizations can likewise be devised in a straightforward fashion for 
25 dominant subspaces with any number of dimensions. 

It should be noted that like the MMSE combiners, the eigenfilter basis itself is also 
implemented using a series of taps in a space-time arrangement. However, two factors 
conspire to make the amount of work involved in adapting these taps negligible. 

Since the eigenfilters can be considered fixed in the short-term, the adaptation takes 
30 place in a very long-term context compared with the MMSE combiner taps. The proposed 
adaptation scheme here is based on the power iteration method. 



■ 
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♦ Only one basis is required to accommodate a plurality (up to A/H) of MMSE 
combiners. 

INDUSTRIAL APPLICABILITY 
5 Embodiments of the invention would be useful in receivers in stations that receive 

multiple signals simultaneously, such as (i) base stations in cellular communications systems 
or access points in wireless LANs; (ii) relay stations or terminal stations in ad hoc or 
unlicensed or packet radio networks capable of maintaining multiple links simultaneously; and 
(iii) terminal or access points of multiple-input multiple-output (MEMO) systems. 

10 Embodiments of the present invention may provide a less costly solution in terms of 

the processing power, the hardware complexity, or both. In fact, they can provide a 
reduction in complexity of an order of magnitude with respect to a canonical linear 
space-time receiver, yet with minimal performance degradation. 

It should be appreciated that the present invention is not limited to the foregoing 

15 embodiments but could be applied equally well in other cases where continuous tracking 

■ 

(using adaptive algorithms such as LMS) is necessary. 

The reduced complexity aspect of preferred embodiments of the present invention 
stems from (i) the shared nature of the common preprocessing section, i.e. it is reused for all 
users; and (ii) the fact that it is adapted slowly, i.e. is less demanding in terms of hardware 
20 and / or software complexity. 

It will be appreciated that the invention is not limited to receivers employing space- 
time processing but embraces receivers employing space-frequency processing, for example 
using Fast Fourier Transforms, or even strictly spatial processing. 

Although embodiments of the invention have been described and illustrated in detail, 
25 it is to be clearly understood that the same are by way of illustration and example only and 
not to be taken by way of the limitation, the scope of the present invention being limited only 
by the appended claims. 

m 

* 

Definitions 
30 In this specification: 

A n channel M refers to the relationship between a transmitted signal and a 
corresponding received signal. 



« 
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A "vector channel" refers either to a channel with a single input and multiple outputs 
(SIMO) or a channel with multiple inputs and a single output (MISO). Each entry in a 
channel vector describes the amplitude and phase of the corresponding channel component. 
A "dispersive channel" or "wideband channel" is a channel with an impulse response 

5 significantly longer than a symbol period, thus resulting in overlap between subsequent 
transmitted symbols, i.e. intersymbol interference (ISI). Such a channel cannot be described 
adequately with a single complex gain. It can either be described as a continuous or discrete 
(i.e. symbol-spaced samples) function of delay in the time domain, or as a continuous 
function of frequency in the spectral domain. 

10 A "narrowband channel" or "flat fading channel" has an impulse response shorter than 

a symbol period and can thus be described by a single complex gain (or a vector of the same 
in the case of a vector channel). 

A "space-time channel" is a vector channel which is also dispersive and can thus be 
described as a space-time matrix of complex gains, or a frequency-dependent vector. 

15 A "covariance matrix" is a time or frequency-averaged outer product of a vector 

quantity or a matrix quantity. In the context of this disclosure, either quantity is either a 
vector signal or a vector channel The covariance matrix of a wireless channel characterizes 
the long-term statistics which underlie the rapid and random instantaneous fluctuations typical 
of such channels. The fact that such covariance matrices vary at a much slower rate than the 

20 channels they characterize is exploited by embodiments of the present invention to reduce 
their complexity. 

A channel "subspace" is a multidimensional space made up of a subset of the 
dimensions making up the AT-dimensional space characterizing an TV-element channel vector. 
A "dominant subspace" is a subspace corresponding to the most significant dimensions of the 
25 channel, i.e. a subset of orthogonal directions which, on average, contain most of the 
channel's energy. 
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